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ABSTRACO? 


Each  service *s  material  probably  will  be  carried  by  rail, 
truck,  ship  or  airplane  while  being  distributed  from  factory  to  the 
ultimate  consumer.  Despite  similarity  of  total  vibration  environ¬ 
ment,  at  least  seven  different  vibration  tests  are  prescribed  by  the 
government,  Hie  two  reasons  for  this  situation  are:  (l)  there  is 
no  basis  for  condensing  long  terra  transport  conditions  into  a  short 
terra  laboratory  test;  and  (2) ,  valid  description  of  Ihe  shipping 
vibration  environment  does  not  exist. 

Emphasis  is  placed  on  developing  a  theoretical  basis  for  a 
laboratory  test  since,  if  for  no  other  reason,  theory  will  define 
the  numerical  values  of  the  shipping  environment  which  should  be 
collected.  Because  a  shipment  may  occupy  a  long  time  period 
during  which  specific  conditions  may  occur  almost  at  random,  use 
of  modern  statistical  theory  (in  this  case  the  theory  of  random 
vibration)  is  indicated.  Equations  are  derived  showing  content 
mean  square  response  in  terras  of  suspension  natural  frequency  and 
damping  ratio  and  the  environmental  characteristics  of  acceleration 
spectral  density  and  bandwidths. 

Linear  fatigue  accumulation  concepts  are  then  used  as  a 
basis  for  realistically  shortening  test  times,  Ihe  literature  not 
being  deer  concerning  whether  standard  fatigue  data  (S-W  data) 
should  be  fitted  to  serai -log  or  to  log-log  curves,  two  expressions 
are  derived  v/hich  equate  sinusoidal  test  response  amplitude  and 
time  id-th  field  time  and  the  field  response  amplitude  previously 
derived.  Hie  equations  show  that  it  is  sufficient,  when  evaluating 
the  suspension,  to  vibrate  at  a  single  frequency,  preferably  the 
lowest  resonant. 

Ob  avoid  recommending  Impractical  tests,  the  capacities  of  exis¬ 
ting  testing  equipment  are  reviewed.  Using  the  results  of  this 
review  as  a  guide,  an  equation  relating  one  sinusoidal  test  with 
another  is  written  using  the  principle  of  equal  input  work.  This  ex 
pression  has  the  same  form  as  the  derived  theoretical  expression 
based  on  plotting  S-N  curves  on  log-log  paper.  Hence,  it  is  not 
essential  to  insure  absolute  con^liance  with  a  stipulated  test  arapli 
tude  provided  test  time  is  varied  accordingly. 


By  use  of  the  foregoing  theoretical  and  practical  approaches 
an  interim  test  procedure  is  recommended  for  standard  use  until 
such  time  as  enough  field  vibration  data  can  be  analyzed  to  provide 
a  better  basis  for  a  standard  test.  Further  recommendations  in¬ 
clude: 


a.  Every  effort  be  made  to  salvage  as  much  random 
data  from  existing  records  as  possible. 

b.  Future  records  be  analyzed  for  contribution  to  a 
proper  statistical  description  of  the  environment. 

c.  Experiments  be  conducted  to  verify  the  presumption 
that  cumulative  fatigue  hypotheses  apply  to  package  cushion¬ 
ing  devices  and  materials. 


INTRODUCTION 


All  progressive  package  designers  agree  that  there  is  no 
way  in  which  the  service  performance  of  a  new  design  can  be 
predicted  in  the  laboratory  without  including  some  simulation  of 
the  vibrations  encountered  in  shipments  This  statement  is  as  true 
of  commercial  packages  as  it  is  of  military  packages.  In  fact  it 
is  standard  practice  to  precede  handling  tests  with  a  vibration  test 
if  onlljr  to  measure  the  performance  of  the  container  in  a  more 
realistic  configuration.  When  so  used,  the  test  .can  be  considered 
as  a  part  of  the  conditioning  process. 

Vibration  isolation  in  missile  shipping  containers  calls 
for  maximum  skill  on  the  part  of  the  container  designer  partic¬ 
ularly  when  one  considers  that  a  suspension  system  also  must  per¬ 
form  as  a  shock  isolator.  Designing  to  provide  for  satisfactory 
isolation  of  large,  infrequently  occurring,  transients  and  smaller 
frequently  occurring  "steady  state"  phenomena  demands  a  compro¬ 
mise  between  design  goals  which  are  essentially  incompatible  except 
in  special  cases. 

While  a  growing  number  of  thoughtful  writers*  believe  that 
primary  emphasis  should  be  placed  on  designing  for  vibration 
isolation  and  letting  shock  isolation  be  the  secondary  goal,  the 
historical  emphasis  in  military  specialty  container  work  has  been 
placed  on  the  shock  isolation  aspects.  Nevertheless,  considerable 
attention  has  been  given  to  the  vibration  control  problem  and  all 
services  have  used  vibration  tests  to  evaluate  shipping  containers. 
The  picture  is  complex,  frequently  confusing  and,  in  many  cases, 
contradictory.  This  report  is,  therefore,  devoted  to  an  analysis 
of  the  problem  and  development  in  detail  of  the  reasons  underlying 
the  recommendations  to  be  made.  By  careful  exposition  of.  the 
reasoning  behind  a  somewhat  arbitrary  recommendation,  we  aim 
at  promoting  intelligent  discussion  and  decision  . 


*See,  for  example,  references  (1)  and  (2).  Throughout  this  report 
numbers  in  parentheses  refer  to  literature  cited  at  the  end. 


DOUBLE  AMPLITUDE  (2y.)  (IN  INCHES) 


FREQUENCY  (IN  CPQ  f 


All  vibration  tests  consist  of  suitable  combinations  of 
frequency,  amplitude  and  time.  In  a  sinusoidal  test,  amplitude 
is  controlled  by  stipulating  linear  excursion  or  by  expressing 
the  number  of  multiples  of  the  acceleration  of  gravity  obtained  by 
relating  acceleration,  frequency  and  linear  amplitude  with  the 
familiar  approximation: 

G  =  ,1025  Cd 

where  G  —  number  of  multiples  of  acceleration  of 
gravity  (386  in/sec.  /sec.  ) 

^  ^  amplitude  (distance  from  rest  position  to 

peak, )  inches 

-P  ^  frequency,  cycles  per  second 

Figure  1  has  been  prepared  for  ready  reference  in  comparing 
these  sometimes  confusing  relations. 

As  of  this  writing,  shipping  container  tests  required  by 
the  government  are  all  of  the  sinusoidal  type.  Study  discloses, 
however,  that  there  are  almost  as  many  tests  as  agencies  capable 
of  imposing  the  test.  Figure  2  compares  the  seven  more  common 
tests.  This  figure  does  not  include  those  previously  required  by 
various  Bureau  of  Ordnance  field  activities  prior  to  issuance  of 
Specification  MIL-W-21927  since  such  a  catalogue  has  appeared 
elsewhere  (3).  Nor  does  the  figure  include  those  tests  which  are 
primarily  designed  to  test  equipment  (e.  g.  MIL.- STD- 1 67,  MIL¬ 
T-18404,  MIL-E-5400,  MIL-T~945,  MIL-STD- 353).  The  pro¬ 
cedures  of  MIL -E- 5272  and  MIL-E-4970  (also  intended  to  be 
equipment  tests)  have,  however,  been  included.  They  sneak  in 
as  shipping  container  tests  in  a  manner  often  unsuspected  by 
those  responsible  for  package  design  and  development. 

We  are  aware  that,  in  the  ballistic  missile  program,  the 
phrase  ’'transit  case"  in  MIL-E-4970  paragraph  4.  6.  8  has  some¬ 
times  been  interpreted  by  Space  Technology  Laboratories  as  meaning 
shipping  container  even  though  such  is  not  sound  packaging  terminol- 
logy  and  in  spite  of  the  fact  that  paragraph  4.  6b  specifically 
eliminates  "equipment  prepared  for  logistic  or  supply  shipment 
and  delivery,  " 
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The  reader  should  also  note  that  Air  Force  Specification 
MIL-P--9024A  vibration  test  information  is  given  in  Figure  2.  The 
latest  version  (MIL-P  - 902 4B)  has  abandoned  aiay  pretense  of  a 
vibration  test  but  relies  on  a  shipping  test  with  a  "suitable"  vibration 
test  as  an  alternate.  Single  sample  shipment  testing  is  time  con¬ 
suming  and,  because  of  the  small  sample,  of  exceedingly  dubious 
validity. 

Now  any  rational  test  could  be  the  basis  for  a  standard  test. 

To  determine  standardization  potential,  one  must  answer  two 
questions: 

(a)  Is  standardization  necessary  or  desirable? 

(b)  Is  standardization  possible ? 

The  answer  to  (a)  is  emphatically  yes.  With  minor  exceptions 
all  missile  containers  are  being  tested  for  exposure  to  essentially 
the  same  environments,  i.  e.  vibratior.s»  encountered  when  shipping 
via  truck,  rail,  ship  or  air.  There  is,  therefore,  no  defensible 
excuse  for  knowingly  continuing  a  situation  which  seems  to  have 
produced  as  many  test  procedures  as  technicians.  In  connection 
with  another  project,  Reed  Research  received  11  replies  concerning 
the  desirability  of  vibration  test  standardization.  Of  these  11,  9 
urged  test  standardization,  1  had  minor  reservations  and  1  indicated 
no  experience. 

The  answer  to  question  (b)  must  be  given  in  two  parts 

{1)  Standa  rdizaticin  now  is  possible  on  an  arbitrary  basis 
if  certain  postulates  be  act  epted. 

(2)  It  appears  pcissible  to  devise  a  standard  form  of  test 
with  quantitative  confidence  that  if  will  predict  service  performance 
provided  certain  further  investigations  are  undertaken. 

The  remairder  of  this  report  is  concerned  with  simultaneous 
development  of  both  of  these  themes  with  the  hope  that  the  test 
suggested  for  use  now  t  an  be  simply  modified  to  become  the  "ideal" 
test  suggested  by  theory. 
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THE  VIBRATION  ENVIRONMENT 


Such  data  as  are  available  indicate  that  the  environment 
includes  the  amplitudes  and  frequencies  shown  in  Figure  3,  This 
figure  is  based  on  the  summaries  contained  in  RD  219/3  (4)  with 
the  significant  exception  that  the  values  given  by  Guins  and  Kell  (5) 
for  AAR  1915  trucks  have  been  omitted.  Reason  for  this  omission 
is  the  very  small  number  of  such  cars  to  be  foxmd  in  interchange 
service  today  and  the  reasonable  presumption  that  any  traffic  man<- 
ager  can  stipulate  not  receiving  such  cars  for  missile  shipments 
even  in  dire  emergency. 

Some  measures  of  statistical  occurrence  have  been  attempted 
as  exemplified  by  Figure  4  taken  from  Ott  (6),  Ott's  curves  are 
cumulative.  That  is,  a  point  on  the  curve  at  any  given  acceleration 
nas  an  ordinate  representing  tne  number  of  peaks  of  that  intensity 
or  greater  which  may  be  expected  duri.ig  the  hypothetical  journey. 

The  existing  data  are  subject  to  certain  devastating 
criticisms: 


a.  Information  concerning  position  of  pickups  and  load¬ 

ing  conditions  of  carriers  is  signally  lacking  in  practically  all  re¬ 
ports. 


b.  Instrumentation  descriptions  are  sketchy  and  state¬ 

ments  concerning  methods  of  data  reduction  leave  much  to  be 
desired. 


c.  The  data  on  probable  number  of  occurrences,  such 
as  summarized  by  Ott,  contain  no  data  concerning  the  frequencies 
associated  with  the  "shocks’'. 

d.  The  data  are  presented  as  samples  reduced,  for 
the  most  part,  as  though  the  inputs  were  complex  waves  and  made 
up  of  discrete  frequencies  and  amplitudes. 

The  last  mentioned  "defect"  is  the  one  that,  apparently, 
places  the  entire  subject  of  shipping  container  vibration  testing, 
and  the  field  data  upon  which  a  test  is  to  be  based,  in  jeopardy. 
Ever  since  the  challenging  paper  by  Morrow  and  Muchmore  (7) 
the  majority  of  vibration  workers  have  been  insisting  that  the 
shipping  environment  is  random  and  that  attempting  to  predict 
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behavior  on  the  basis  of  classical  vibration  theory  simply  does 
not  correlate  with  experience. 

Now  classical  theory,  in  effect,  tells  us  that  if  we  can 
identify  all  the  components  of  a  vibration  record  in  terms  of 
certain  relatively  simple  functions  (sinusoidal,  step,  etb.  )  we 
could  predict  behavior.  In  the  case  of  transportation  vibrations, 
however,  there  are  so  ma©y  variables  affecting  the  vibration 
that  it  seems  to  many  thoughtful  people  that  it  is  impossible  to 
predict  precisely  what  will  happen  in  the  next  instant  of  time’!'. 
Probability  theory  can  be  used,  however,  to  predict  what  might 
reasonably  be  expected  to  happen  over  a  given  period  of  time  on 
the  assumption  that  any  of  the  possibilities  could  happen  next. 

In  effect,  we  treat  the  problem  much  as  we  would  treat  the 
problem  of  experimental  error.  Hence  the  phrase  "random 
vibration". 

Mathematicians  call  the  models  derived  to  handle  such 
problems  stochastic  processes**!'.  The  stochastic  process  is 
used  in  much  the  same  way  as  a'simple  sinusoidal  description 
is  used.  It  must  be  remembered  that  both  are  simplifications 
of  reality  but  the  stochastic  is  considered  a  better  simplification. 

In  fact,  Morrow  (11)  states  "For  mis*siles  the  most  common  and 
most  damaging  tactical  vibration  is  random". 

As  implied,  the  concept  of  random  vibration  involves  appli¬ 
cations  of  mathematical  statistics  to  mechanical  problems.  This 
is  not  altogether  a  new  concept.  Einstein's  (12)  1905  theory  of 
Brownian  motion  (which  may  be  considered  a  special  case  of 
random  vibration)  was  extended  to  strings  and  certain  beams  in 
1931  (13).  In  1933  Wiener  (14)  introduced  such  concepts  as  power 
spectral  density  (to  be  defined  later)  in  his  generalized  harmonic 
analysis  giving  a  powerful  analytical  tool.  The  big  impetus  towards 
current  use  of  random  theory  was  given  by  S.  O.  Rice  in  his  classic 
papers  on  random  noise  in  electric  circuits  (15),  (16).  A  complete 
summary  of  available  knowledge  is  beyond  the  scope  of  this  work. 

Of  particular  interest,  however,  are  the  analyses  of  random 


*For  much  of  what  follows  concerning  history  and  definitions  we 
are  indebted  to  the  clear  exposition  by  Crandall  (8,9)  and  by  James 
Nichols  and  Phillips  (10). 

**Stochastic  could,  if  one  wished,  be  considered  a  ten-dollar  word 
implying  chance. 
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FIGURE  4 


DIAGRAMS  SHOWING  INTENSITY - 
STATISTICAL  FREQUENCY  OF 
VIBRATION  IN  TRANSPORTATION 
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vibrations  in  beams  (13),  (17),  (18),  (19),  in  plates  (17),  (20) 
and  in  shells  (21). 


The  statistical  approach  to  vibration  analysis  requires  the 
introduction  of  certain  concepts  which  may  be  new  to  those  grounded 
in  classical  vibration  concepts.  For  clarity,  we  reintroduce  these  con- 
copts,  basing  our  exposition  largely  on  references  (9)  and  (10). 


Let  us  consider  any  motion,  %  (•t')  i  in  which  y.  denotes  a 
displacement  and  b  indicates  that  X  is  a  fimction  of  time.  This 
motion  may  be  the  motion  of  the  bed  of  a  truck  or  railroad  car  or 
airplane  or  other  support  for  a  package  and  represents  the  forcing 
function  in  the  dynamical  problem  under  study.  We  will  designate 
by  ^  the  absolute  acceleration  of  this  motion. 


Three  statistical  measures  are  of  particular  interest  with 
respect  to  this  motion.  The  first  is  the  mean  square  defined  by 

which,  as  the  name  implies,  is  a  very  general  measure  of  intensity. 

Sometimes  the  square  root  of  the  mean  square,  referred  to  as  the 

root  mean  square  (and  also  as  the  standard  deviation),  is  used  and 

this  is  ^  ^ 

^  -  g  A 


j 


These  average  inXensity  measures  give  no  indication  of  the 
oscillatory  nature  of  the  motion.  This  is  pictured  by  the  power 
spectral  density,  S  which,  as  the  name  implies,  shows 

the  power  content  of  the  vibratory  components  as  a  function  of 
their  circular  frequencies,  OO  - 

The  power  spectral  density  has  the  property  that  its 
integral  is  the  mean  square  of  the  motion.  That  is 

X  =  ihyr  f  S  (ro^  d  CO  BJ 

This  simply  means  that  the  average  power  in  the  motion  is  the  sum 
of  the  contributions  of  each  infinitesimal  component  at  each  frequency. 

It  is  considered  characteristic  of  random  motions  that  power  spectral 
density  is  a  boiinded  function.  In  the  case  of  transportation  vibration 
it  is  also  considered  reasonable  that  the  function  will  contain  peaks 
in  the  neighborhood  of  the  resonant  frequencies  of  the  vehicle. 
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The  power  spectral  density  is  also  defined  mathematically 
in  terms  of  the  Fourier  Transform.  Thus,  *if  we  consider  a  function 


X  ct)*xct)  when  -  T<ii  <  T  with  T  a  particular  value  of 


time  and  X  (t)=  0  otherwise 


and  write  A  CO) - 

T 


X  Cc)  e 

•T 


dt  * 


-i,C0T  ,, 

X  ct;e  dt 

T 


-T 


then  JP is  the  Fourier  Transform  of  .  If  X(t)  is 

considered  to  extend  indefinitely  on  the  time  scale,  then  this  limit 
does  not  exist.  However,  the  limit  jLlifTl  -L  I A  exist,  if 

T-*  «>o  ^ 

X(t)  is  a  random  motion,  and  this  limit  defines  the  power  spectral 
density.  That  is 


At  (a) 


? 


Such  power  spectral  density  functions,  when  plotted  against 
the  frequency,  often  show  marked  peaks  at  certain  frequencies,  A 
useful  description  of  these  peaks  is  given  by  the  peak  value,  the 
frequency  at  the  peak  and  the  band  width.  The  latter  is  conventionally 
measured  by  the  frequency  difference  between  the  points  on  each 
side  where  the  spectral  density  has  fallen  to  one  half  the  peak  value. 


The  third  important  statistical  measure  of  random  vibration 
is  the  probability  density  function.  The  underlined  phrase  is  drawn, 
at  least  in  part,  from  the  language  of  mechanical  engineering.  Be~ 
cause  of  mathematical  similarity,  it  is  convenient  to  consider  that 
the  probability  of  a  value  falling  within  a  certain  bounded  region  is 
equal  to  the  mass  of  that  region.  It  is  usual  to  assume  that  the 
probability  density  of  a  random  motion  is  normal  or  Gaussian.  That 
is  to  say,  the  probability  of  encountering  a  particular  value  of 
is  as  shown  in  Figure  5. 


’i'Alternati vely,  consider  the  autocorrelation  function 

R^(T)  ir  * 

Then  the  power  spectral  densTty  is 

f  “  r_^(t)  dfc 

ryur* 

That  is,  the  power  spectral  density  is  twice  the  Fourier  Transform 
of  the  autocorrelation  function. 


II 
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"When  the  response  of  a  resonant  system  to  a  random  exci¬ 
tation  is  considered,  an  envelope  of  peak  values  can  be  drawn.  If 
the  excitation  is  Gaussian,  and  the  a  of  the  system  is  high*,  then 
the  probability  density  function  for  the  envelope  is  a  Rayleigh 
density  as  shown  in  Figure  6.  -Xp 


The  equation  of  this  curve  is: 


P(*p) 


(6) 


where; 

p(xp)  is  the  probability  density  of  xp;  i.  e.  ,  is 

the  probability  thaty^  lies  between  Kp  andXpVAxp 

Xp  is  a  peak  value  of  X  or  an  ordinate  of  the  envelope 
curve, 

Cy  is  the  rms  value  of  X  . 

It  is  clear  frona  the  foregoing  that  an  adequate  description  of 
the  environment  useful  in  a  statistical  approach  requires  data  con¬ 
cerning  the  envelope  of  all  possible  acceleration  spectra  together 
witfi  a  measure  of  the  bandwidths  of  the  peaks  therein  and  an  approxi¬ 
mate  measure  of  possible  duration.  It  is  also  clear  that  currently 
available  information  concerning  vibrations  encountered  in  shipment 
cannot  be  considered  a  statistical  description  of  the  environment. 

^  About  all  that  can  be  said  about  these  data  is  that  the  recorded  fre¬ 
quencies  do  occur  and  that  the  reported  amplitudes  occur  a  sufficient 
number  of  times  to  be  noticeable.  Statistical  measures  are,  however, 
still  insufficient  for  completely  logical  test  design.  In  effect,  therefore, 
no  matter  what  test  is  suggested  at  present  writing,  it  must  be 
accepted  that  no  quantitative  confidence  figure  can  be  derived  con-  , 
cerning  faithful  prediction  of  service  performance. 

The  impossibility  of  deriving  a  quantitative  confidence  level 
should  not  deter  the  government  from  selecting  a  ^'standard*'  test 


*That  is,  to  say,  the  system  is  lightly  damped.  By  definition 


^  1  Cc 

Co  where  J 


damping  ratio, 


CjjS*  damping  coefficient,  lb,  sec.  /in, critical  damping  coefficient, 
«  ssspring  rate  of  suspension  and  TTl  mass. 
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based  on  qualitative  confidence  levels  derived  from  experience 
provided  it  is  recognized  that  an  arbitrary  "standard”  has  been 
established  and  that  work  should  continue  towards  achieving  a 
better  standard.  The  major  difficulty  with  selecting  a  standard 
now  is  that,  objectively,  each  method  now  in  use  has,  with  certain 
exceptions,  equal  standing  with  every  other  method,  including 
the  possible  standard  to  be  suggested  here. 

Our  problem  is  to  establish  a  theoretical  basis  for  a 
transportation  vibration  test  assuming  that  a  valid  statistical 
description  of  the  environment  exists  and  then  examining  what 
practical  compromises  are  necessary  to  permit  performance  of 
any  test  now.  Even  to  go  this  far  we  must  proceed  from  two 
assumptions  which  we  will,  hereafter,  regard  as  axiomatic: 

a.  Shipping  container  vibration  tests  evaluate  only 
shipping  containers.  In  other  words,  one  is  testing  performance 
of  the  isolator  system,  not  performance  of  the  missile.  "Trick" 
frequencies  and  amplitudes  are  a  problem  for  the  missile  designer. 
If  they  exist  in  the  range  to  be  considered,  supplementary  require¬ 
ments  may  be  written  but  not  by  packaging  personnel. 

b.  Any  test  is  a  simulation  not  a  duplication  of  the 
environment.  In.  a  laboratory  test,  the  time  scale  is  reduced 
and,  most  frequently,  severity  is  increased  to  compensate. 
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MODEL  OF  A  PACKAGED  ITEM 


RESPONSE  OF  A  PACKAGED  ITEM 
TO  RANDOM  MOTION 


With  the  foregoing  assumptions  in  hand,  we  proceed  to  an 
analysis  of  the  response  of  a  package  with  a  linear  damped  spring 
system  such  as  shown  in  Figure  7.  For  simplicity  we  will  restrict 
ourselves  to  the  linear  case.  *  If  the  system  be  linear  then,  by 
use  of  the  concept  of  normal  modes,  it  is  possible  to  treat  the 
motional  response  of.  the  item  with  respect  to  the  container  as  a 
single  degree  of  freedom  system. 


The  input  motion  "coordinate"  is  ct)  with  corresponding 


.e 

acceleration  Q.  The  absolute  response  coordinate  is 

Xo  (t)  and  the  relative  motion  The  absolute  response 

acceleration  is  Q  (jj)  ^  It  must  be  recognized  that  the 


absolute  response  may  require  a  vector  pair  to  describe  its 
translational  and  rotational  components.  The  response  coordinate 

XoCt)  is  simply  a  magnitude  which  characterizes  the  amplitude 


of  the  response  in  the  normal  mode  under  consideration-  The  same 
thing  is  true  of  the  relative  motion*  So  far  as  the  input  is  concerned, 
it  can  generally  be  assumed  to  consist  of  a  rectilinear  motion,  although 
in  principle  it  could  just  as  well  be  a  rotation  or  a  combined  motion. 

If  X  ,  is  a  complex  exponential  function  of  time,X  ' 
then  the  output  is  ^ 

and  the  relative  motion  is  X  = 


*For  delicate  items  the  mounting  most  often  consists  of  shear  mounts 
which  are  substantially  linear  in  the  range  under  consideration.  This 
restriction  does  not,  therefore,  result  in  an  academic  enterprise. 
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,  are  complex 


where,  in  general, X,  ,X  and  X- 

u  ’  o 

functions  of  the  circular  frequency,  00 


The  ratios  (22) 


X. 


X. 


and 


X. 

V 


(«-y 


.in  which  H  is  the  complex  frequency  response  or  magnification 
factor,  permit  us  to  write 


and 


o  o  , 


where  is  the  natural  frequency  of  the  normal  mode  under 
consideration 


(7) 


(8) 


(9) 

(10) 


and 


is  the  ratio  of  the  damping  to  the  critical  damping. 


When  the  input  acceleration,  01 .  ,  is  given  the  response 
acceleration  is  O,^  and  ** 

Qq  ^  Xo  „T  (11) 

a.  "■  X.  '  ' 

U  U 

The  relative  motion,  ,  is  given  by 


From  these  relations  between  excitation  and  the  response 
thereto  it  is  possible  to  obta^  relations  between  the  statistical 
■measures  of  each.  Thus, if  ^  ^^^(co^enote  the  power  spectral 


density  of  the  input  motion  and^^^(^CAi^nd  on  X  the  spectral 

densities  of  the  absolute  and  relative  responses,  respectively,  then  (23) 


Sox  =  I  Tq!  Si  (.Co)-  1-I-4T  _ Six  (w) 
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(14) 


S  (co)  = 


n 


S.  ("'l 

LX 


V 

Leo  h .) 


[K- 


As  the  input  is  more  apt  to  be  given  in  terms  of  acceleration  and 
as  the  absolute  acceleration  response  is  of  more  interest  than  the 
absolute  response  motion,  we  denote  the  input  acceleration  spectral 
density  by  S>viQ,  (CO)  and  the  acceleration  spectral  density  of  the 
absolute  response  (CO^, 

Then 

Soc,(''-^'^=lTI^S.  CcS) 

S  fo) 

y  h  vj  cx^ 


(15) 

(16) 


From  these  the  mean  square  absolute  acceleration  response  can 
be  determined  as: 


_ a. . 

o.„ 


1  X  \"'SiJiS)dc^ 


and  the  mean  square  relative  motion  response  as: 

U 

a*  '  1  H 


<3  u  OJ 


If  the  input  acceleration  spectral  density  is  flat  in  the 
neighborhood  of  the  natural  frequency,  then,  approximately. 


and 


Q.0 


X, 


.3. 


1 


Where  ^  0L.CCOj^)ienotes  the  value  of  the  input  acceleration 


(17) 


(18) 


(19) 


(20) 


density  when  the  input  frequency  is  equal  to  the  natural  frequency 
of  the  suspended  mass.  If  the  input  acceleration  spectral  density 
has  a  peak  at  CO  -  CO  ir,  ,  as  is  to  be  expected,  Soq_  (jX)^  is  a  variable 
during  the  integration  of  equations  17  and  18  .  The  following 

function  is  considered  descriptive  of  the  variation  of  : 


This  function  requires  for  its  specification  only  its  peak  value, 

the  corresponding  frequency,  CO  n  ,  and  its  bandwidth, 

coh  6c.  which  is  the  frequency  difference  between  the  points 

on  each  side  of  the  peak  at  which  the  ordinate  of  the  curve  is  half 
the  peak  value.  Putting  equation  into  equation  and  inte- 
grating  {see  Appendix  I)  gives  the  means  square  of  the  output 
acceleration  in  terms  of  the  natural  frequency  of  the  suspension 
system  and  its  damping  ratio,  as  well  as  the  input  acceleration 
density  at  CO  3  6D  ^  and  its  bandwidth: 


Similar  substitution  in  equation  ps]  produces  the  following  expres 
Sion  of  the  mean  square  relative  motion  response: 


SINUSOIDAL  STEADY  STATE 
versus 

RANDOM  VIBRATION  TESTING 


Equipment  for  producing  random  motion  in  the  laboratory 
is  very  expensive  and  not  readily  available  to  all  who  might  be 
involved  in  the  packaging  of  complex  items.  On  the  other  hand 
sinusoidal  type  testing  machines  and  their  associated  devices  and 
controls  are  considerably  less  expensive.  If  they  could  be  used  a 
meaningful  test  would  be  within  the  reach  of  a  greater  number  of 
operators  and  testing  costs  {hence  the  cost  of  packaging)  could  be 
greatly  reduced.  We  now  address  ourselves  to  this  problem. 

If  a  sinusoidal  test  motion  with  amplitude  is  applied  at 

the  natural  frequency  of  the  package,  a  uniform  re spons ^amplitude 
will  result.  The  mean  square  value  of  this  response,  n/  , 

is  given  by  \/  ^  . 


(24) 


For  one  form  of  equivalence,  mean  square  responses  in  the  test  and 
in  transportation  could  be  made  equal.  Combining  [24]  and  |2^ 


)<;t 


(it  "If) 


V  B  i,  '  A 

The  test  acceleration  amplitude^ for  such  equivalent 
motion  will,  since  be  given  by 


B;, 

Such  a  test  would  involve  power  dissipation  and  heat 
generation  within  the  package  substantially  equal  to  that  in  service. 
It  would  fail  to  simulate  service  conditions,  however,  with  respect 
to  the  distribution  of  amplitudes.  For  example.  Figure  8  shows 
that,  if  the  distribution  of  amplitudes  in  service  is  a  R.ayleigh 
distribution,  the  sinusoidal  test  amplitude  (which  is  times  the 
rms  value)  is  exceeded  40%  of  the  time  in  service.  A  service 
amplitude  twice  the  sinusoidal  teat  amplitude  is  exceeded  about 
1%  of  the  time  and  thrice  the  test  amplitude  0.  01%  of  the  time. 


(25) 


(26) 
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FIGURE  8.  PROBABILITY  OF  EXCEEDING  ANY  VALUE 
OF  PEAK  DISPLACEMENT,  IF  THE  PEAKS 
FOLLOW  THE  RAYLEIGH  PROBABILITY  DENSITY 
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The  possibility  exists,  therefore,  and  particularly  where  heat  is 
not  the  cause  of  failure,  that  service  amplitudes  in  excess  of  rms 
value  of  the  test  amplitudes  would  occur  a  sufficient  number  of 
times  to  cause  a  fatigue  failure  of  the  suspension. 


It  should  also  be  noted  that  equations  1^253  and  have  no 
test  duration  factors  and  hence  imply  that  the  test  time  must 
approximate  the  transit  time.  ‘Thus,  such  a  method  of  stipulating 
the  test  parameters  fails  on  two  counts.  In  the  next  section,  there¬ 
fore,  we  explore  the  nature  ol  fatigue  failures  with  a  view  to 
determining  what  an  equivalent  test  time  could  be  with  increased 

'^C^or^  above  the  values  indicated  by  equations  [zs]]  and  ^26^. 
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FATIGUE  AND  SHORTENED 
TEST  TIMES 


Shortening  test  time  by  increasing  severity  is  a  procedure 
which  has  been  followed  intuitively  in  almost  all  disciplines  for 
years.  Modern  philosophical  justification  for  so  doing  in  a  vibra¬ 
tion  test  is  based  on  cumulative  fatigue  damage  concepts  developed 
by  Palmgren  (24),  Miner  (25),  Corten  and  Dolan  (26)  and  others*. 
These  hypotheses,  based  on  metal  fatigue,  postulate  that  a  stress 
at  some  level  reduces  the  life  of  the  test  object  by  a  stipulated 
fraction  of  its  fatigue  life  at  that  stress.  Repeating  the  shocks 
produces  cumulative  damage  eventually  resulting  in  fatigue  failure. 
Miner's  original  hypothesis  showed  linear  damage  accumulation 
characterized  by  the  relation 

y  iii-=  C 

K.  ^ 

in  which  n  (^represents  the  actual  number  of  cycles  at  a  given  stress 

level, A/  the  fatigue  life  at  that  same  stress  level  and  ^  is  a. 

constant.  Miner  proposed  that  •  Obut  limited  experimental 
results  show  a  range  for  ^  from  0»3  to  about  3i0,  depending  upon 
whether  the  high-stress  cycles  come  early  or  late  in  the  test* 

Corten  and  Dolan's  non-linear  accumulation  proposal  correlates 
with  experimental  data  on  cold-drawn  wire  and  is  characterized 
by  the  equation  .  / 


where  Iv^represents  the  total  number  of  cycles  to  failure;  Vi  ,  the 
fatigue  life  at  the  maximum  stres  are  the  per¬ 
centages  of  cycles  at  stresses  S|  . S  0  ’  ^ 

experimentally  determined  exponent  which,  for  a  specific  cold- 
drawn  wire,  was  found  to  be  6.  57.  Figure  9  compares  the  progress 
of  damage  at  a  given  stress  level  under  the  two  theories. 


Miles  (28)  advances  the  argument  that  presently  available 
fatigue  data  on  metals  scarcely  warrant  non-linearity  complications 


*McClintock  (27)  gives  a  good  resume  of  niany  of  the  theoretical 
investigations. 
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in  predicting  fatigue  life  of  a  single  sample*  He  further  demon¬ 
strates  that  using  non-linear  approaches  has  no  appreciable  effect 
when  the  stress  amplitudes  are  continually  distributed  over  a  wide 
range.  Mains  (29),  using  Miner's  approach  with^O.^and  Corten- 
Dolan's  withd  =  iJj-Slencountered  the  "interesting  coincidence"  that 
life  prediction  in  one  case  was  substantially  the  same*  Since  the 
transportation  input  vibration  is  wide  ranged  and  of  long  duration 
Miles’  argument  may  well  be  applicable  to  the  transport  condition. 
The  test  condition  is  assumed  at  a  constant  stress  amplitude.  We 
feel  free,  therefore,  to  use  Miner's  hypothesis  as  a  point  of  de- 
parture. 

Miner's  hypothesis  can  be  restated  as 

p  _  V 

where  3  —  stress  level,  which  must  be  greater  than  the 

endurance  limit,  stress  at  wliich  fatigue  life  is  infinite. 

In  an  extended  duration  random  motion  increasingly  finer 
stress  gradations  occur  and  we  are  led  to 

CT^  r Ad. 

J  hN/(^ 

in  which  is  now  considered  the  number  of  stress  cycles'  which 
do  not  exceed  S  . 


Sincelnts  a  function  of  time  and  stress  is  proportional 
to  relative  motion  we  may  write  the  following  expression  for  life 
expenditure  as  a  function  of  time 


L  (Q  * 


.  h 


in  which  is  tlie  number  of  cycles  to  failure  at  an  amplitude^p  , 
It  remains  to  express  In  <And  IVyas  functions  of  V  p  in  order  to 


carry  out  the  integration. 


(31) 


We  have  shown  that  the  relative  motion  response,  r\  , 
is  reasonably  narrow  banded  due  tovthe  peaked  character  of  the 
excitation  spectral  density  curve.  Thus  the  number  of  peaks  occur- 
ring  per  second  is  approximated  by  If  we  further  assume 

that  the  suspension  system  has  a  high  Q  ,  i.  e.Q^l  0,  then  the 
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number  of  peaks  is  n  per  second.  *  The  number  of  such 

occurrences  for  which  the  peak  value  lies  between^  p  and  -t-ciy  p 
is  (30): 


—  y 

r|  D  -  ^  ^  J  V  n 


^  —-Si* 

in  which  response  motion  {  ®ri"^  n.  5 

'^ps'peak  value  of  y^ 

The  form  of  the  relation  between  maximum  stress  in  a 
fatigue  test  and  number  of  cycles  to  failure  is  expressible  as  the 
semi-logarithmic  function  ^  ^  ^ 


(32) 


S  r  oC  -  /?  hx^ 


N 


S>5e  (33) 


which  in  terms  of  relative  motion  in  the  package  can  be  written  as 

Vp  =  4  “  0  fv  / 


or 


-^  =  e 

IV. 


u 

-i  ^ 

8  e  s 


X  p  >  - 


(34) 


where  A  and  0  ^re  empirical  constants  (  (^corresponding 
to  the  slope  of  the  3“  curve)  andX  <o  is  the  amplitude  corresponding 

to  Se 


Substituting  ^  32*^  and  Q4]into  ^3^  and  performing  the  indicated 
integration  (see  Appendix  1)  results  m 


9*  ^  A 

£  (t)  =  .^ilL  e 


a 


£JL 


L.  6  J 


a 


-  CP 


fwV' 


^Since  packages  usually  have  a  a<  5  this  appears  a  sweeping 
assumption  but  it  is  not  necessarily  crippling  to  what  follows.  Failure 
can  be  intellectually  divided  into  heat  failure  and  fatigue.  If  heat  is 
the  cause  the  temperature  will  build  up  rather  rapidly  until  some 
equilibrium  is  reached.  This  will  or  will  not  cause  failure-  If  the 
system  has  a  very  low  undamped  natural  frequency,  auxiliary  dampers 
will  be  used.  These  are  expected  to  fail  first  and  then  we  do  have  a 
high  Q  package. 


in  which 


On.  — ^^^value  of  relative  motion  response  to  the  field 


environment  defined  by  equation 


1/  -  Xc 


GrF (l‘^a;he  probability  integral  J  S. 

■^0 

which  can  be  obtained  from  tables,  (31) 


values  of 


Miles  (28)  has  suggested  that  ".it  seems  reasonable  to 

ignore  the  possibility  of  an  endurance  limit  when . stress 

amplitudes  are  distributed  over  a  wi^e  range".  If  this  be  true  in 
fact,  and  it  is  easy  to  visualize  many  distributed  cushioning  materials 
in  which  an  endurance  Iknjt  would  be  very  poorly  defined  or  even  non¬ 
existent,  then  equation  I  35  |  simplifies  to 

E  Ct")  =  e 


(36) 


The  equation  cinalogous  to  equation  [S^fara  sinusoidal 
vibration  test  is  ,\j- 

E.(t)  e(V-t) 


(37) 


where  t  T  =  test  time  and  Xn  “I*'*-  test  response  amplitude. 


To  account  for  possible  variation  from  Miner’s  failure  at  unity 
hypothesis  we  set 

E  .p  (t)  -  c  [  it) 

where  ^  =  the  actual  value  assigned  to  Miner's  constant. 


(38) 


Taking  [35],  ^3'^and^^and  solving  for  time  produc 
O'h.  "•  ^  Q  y  n  T 


Ct 


r_  =  e 


es 


B  '' 


-  Ch 


+  <=■ 


-K,- 


(39) 


or,  substituting  witl^^n  lieu  of 

crZ-XB/nT 


[•'] 


ct 


I-  =  c 


I  4-  ^ 
G 


(40) 
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Miles  (28)  has  also  suggested  that  it  is  equally  convenient 
to  express  fatigue  life  as  a  log-log  relation;  log  S  * 

^s-  If  we  follow  this  through  (see  Appendix  1)  we  arrive  at  the 
following  complete  statement  of  the  ratio  of  test  time  to  service 
life 


Ci,  ■  ^xT  ^  r  (I  +  3^ ) 


(41) 


where  P0+  complete  Gamma  function  which  can  be 

evaluated  from  tables  or  curves  (32). 

incomplete  Gamma  function.  Curves 

.V"  .  . 

are  given  in 


for  some  values  of  the  ratio 


r 


Jahnke-Emde  (32)  while  extensive  tables  are 
given  by  Pearson  (33). 


If  K  be  sufficiently  small  to  permit  ignoring  its  contribution 
41 1  becomes 


-ir  ®  • 


t-t 


'll 


(42) 


It  is  interesting  to  compare  the  effects  of  varying  test 
amplitudes  on  test  times.  For  a  specific  package  design,  the 
extensive  modifiers  will  remain  constants.  Thus,  the  semilog 
equations  can  be  rearranged  to  produce 


T 


while  the  log-log  S-A'^ curve  equations  lead  to 

JlL.=  '/b  . 


(44) 


It  is  appropriate  to  review,  at  this  point,  what  the  theo¬ 
retical  approach  does  and  does  not  tell  us.  The  equations  developed 
herein  indicate  that; 
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a.  For  a  specific  package  a  quantitative  relationship 

between  sinusoidal  test  time  and  service  life  in  a  random  environ¬ 
ment  can  be  computed  if; 

I.  Suspension  systems  follow  the  cumulative 
damage  hypotheses  and  provided,  further,  that; 

(i)  The  slope  of  and  at  least  one  point 

on  the  S.  N^curve  for  various  materials  can  be 
found  experimentally. 

(ii)  The  level  of  the  endurance  limit, 
if  it  exists,  can  also  be  found, 

(iii)  Suitable  values  to  use  for  Miner's 
constant  can  be  found. 

Z.  The  response  characteristics  of  the 
particular  package  can  be  determined.  Note  thatOy^^is 
a  function  of  system  damping  and  input  characteristics, 

3,  Adequate  descriptions  of  the  environment 
are  available.  These  descriptions  must  include  acceleration 
spectral  density,  bandwidth  and  time.  Unless  we  have 
these,  we  cannot  compute (7,^ nor  can  we  compute  1^. 

Since  we  do  not  now  have  the  critical  information  shown 
above,  it  would  aprpear,  at  first  glance,  that  further  progress 
towards  standardization  is  impractical.  This  is  true,  however, 
only  if  we  insist  that  we  must  have,  now,  a  quantitative  service 
life  prediction.  There  are  certain  other  teachings  ir  olicit  in 
equations  m  through  These  are: 

1.  A  sinusoidal  test  can  be  an  adequate 
simulation  of  a  random  environment. 

2.  The  time  of  test  can  be  varied  provided 

there  be  a  concomitant  variation  in  test  response  amplitude. 
The  minimum  such  amplitude  should  be  such  as  to  produce 
a  stress  greater  than  the  endurance  limit.  The  maximum 
is  the  one  time  failure  stress.  In  practice,  other  con¬ 
siderations  such  as  available  clearance  should  probably 
govern. 
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3.  Since  test  severity  is  directly  related  to 
response  amplitude,  it  is  only  necessary,  to  evaluate 
life  of  the  suspension,  to  vibrate  where  response  ampli¬ 
tude  is  the  greatest,  i.  e,  at  the  lowest  resonant  frequency. 

These  three  major  factors  permit  logical  simplification 
of  the  confusion  engendered  by  the  conflicts  inherent  in  Figure 
2.  However,  it  must  be  recognized  that  actual  standardization 
is  only  possible  if  certain  arbitrary  postulates  be  accepted.  In 
addition  to  the  two  postulates  stated  at  the  end  of  the  section  on 
environment  we  will  proceed  on  the  assumptions  that  the  follow¬ 
ing  are  self-evident: 

a.  A  test  written  as  a  specification  requirement  today 
must  be  capable  of  performance  on  equipment  available  today.  If 
not  capable  of  performance  it  dies  from  non-use.  Besides,  how 
could  it  be  a  quasi-legal  basis  for  acceptance  or  rejection  under 
contract? 

b.  Since  it  is  clear  that  some  arbitrary  decisions  must 
be  made,  it  is  possible  for  all  concerned  to  be  arbitrary  in  the  same 
way. 


With  these  thoughts  in  mind,  we  now  proceed  to  a  discussion 
of  the  arbitrary  decisions  to  be  made,  lacing  our  argument  with 
the  lessons  of  the  theoretical  discussion  just  concluded  and  various 
practical  matters  inherent  in  specialty  containers. 
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SELECTION  OF  FREQUENCY  RANGE 
TO  BE  COVERED 

The  theoretical  discussion  has  shown  th^t  we  need  only 
test  at  the  lowest  resonant  frequency.  We  now  examine  which 
particular  frequency  limits  can  be  placed  on  the  specification 
requirements.  There  appears  to  be  no  question  but  that  the 
tests  should  include  the  very^low  frequency  range  of  2-10  cycles 
per  second*  From  the  evidence  available  the  peak  amplitudes 
of  the  random  environment  can  be  expected  in  this  range* 

Further,  container  shock  isolation  systems  for  very  delicate 
components  quite  often  exhibit  natural  frequencies  below  10  cps. 
Such  mounting  systems  are  quite  "soft”,  frequently  require 
supplementary  dampers,  and  resonance  can  produce  large 
magnifications  and  considerable  heat  built-up  which,  in  itself, 
can  damage  the  materials  used  in  mounts  or  cushions. 

Now  RD  219/3  shows  that  higher  frequencies,  up  to  7,  000 
cps,  are  within  the  realm  of  possibility.  The  sinusoidal  amplitude, 
however,  even  at  5g  (an  amplitude  not  claimed  for  the  7,  000  cps 
frequency)  is  extremely  small,  i.  e.  ,  0.  00000102  inches,  which 
is  hardly  worth  considering  to  be  contributory  to  any  response 
motion  of  the  item  in  its  shock  isolation  system.  At  I  00  cycles 
per  second,  amplitude  is  only  .  005  inches  for  5g  accelerations. 

At  60  cycles  per  second,  1*  3g  produces  an  amplitude  of  only  .  003 
inches.  It  is  highly  probable  that  such  amplitudes  would  be  damped 
out  by  almost  any  simple  structure  such  as  a  single  layer  of 
corrugated  fibreboard.  In  a  missile  container,  of  course,  there 
is  considerably  more  structure  than  a  single  layer  of  fibreboard 
between  the  source  of  the  excitation  and  the  suspended  mass.  We 
have  never  personally  favored  higher  frequency  vibration  testing 
of  shipping  containers.*  Because  of  possible  prejudice,  therefore, 
an  industry  survey  went  into  the  question  of  what  is  learned  when 
testing  containers  over  60  cps.  Of  eleven  replies  received  4 


*Some  will  recall  an  early  version  of  MIL-P-7936  in  its  uncoordi¬ 
nated  form  which  was  prepared  in  a  considerably  greater  state  of 
ignorance  by  one  of  the  authors  some  six  years  ago.  Here  testing 
over  55  cps  was  not  required  except  for  rigidly  mounted  components 
in  all-metal  containers. 
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flatly  stated  that  they  had  never  learned  anything,  2  stated  they 
had  no  experience,  2  didn't  answer  this  question.  One  pled  lack 
of  experience. but  doubted  necessity.  One  said^"not  a  great  deal 
although  some  secondary  resonances  were  found  and  a  single 
case  of  poor  attenuation  was  noted.  "  The  last  responder  stated 
that  about  all  he  had  proved  is  that  "you  end  up  testing  the 
resonant  frequency  of  the  fixtures  that  are  supposed  to  transmit 
the  vibration  from  the  table  or  head  to  the  box". 

4 

The  last  remark,  delivered  in  half  joking  language,  truly 
represents  a  basic  facet  of  the  problem.  With  very  few  exceptions 
(the  giant  machine  developed  by  Norair  Division  of  Northrop  or  the 
special  devices,  only  one  of  which  is  currently  known  to  be  in 
existence,  developed  by  Royal  Jet)  there  are  no  machines  which 
will  shake  the  whole  container  at  these  frequencies  if  the  container 
be  much  larger  than  a  case  of  canned  goods.  Some  "cobbled  rig" 
is  almost  always  necessary  to  transmit  force  from  the  relatively 
small  geometry  of  the  shaker  unit's  moving  element  to  the  very 
much  larger  container.  The  higher  frequency  low  amplitude 
vibrations  simply  get  lost  in  this  structure. 

Another  way  of  looking  at  the  utility  of  testing  above  60  cps 
is  to  consider  G-factor  vs.  height  of  drop  which  would  produce  a 
natural  frequency  at  60  cps.  Inasmuch  as,  for  a  linear  system, 

=  (45) 

and 

w  W6) 

°  k.  ^ 

=  .45-5f„JTT-|  (47) 

Figure  10  is  a  plot  of  this  last  equation  for  various  values 
of  'Th  •  This  graph  shows  that,  to  achieve  a  natural  frequency  of  60 
cps  at  a  reasonable  drop  height,  the  contents  must  be  very  rugged 
indeed.  It  is  difficult,  except  in  unusual  circumstances^  to 
envision  stress  levels  on  the  order  of  I /50th  to  l/l60th  of  a 
conservatively  estimated  breaking  stress  resulting  in  fatigue 
damage  to  the  isolator  system. 


K 

therefore 

0 
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RELATION  OF  G“  FACTOR,  HEIGHT  OF  DROP 
AND 

TRANSLATIONAL  UNCOUPLED 
NATURAL  FREQUENCY 
FIGURE  to 


In  other  words,  so  long  as  design  emphasis  is  on  shock 
isolation,  high  frequency  testing  does  not  appear  to  evaluate  any 
significant  life  property  of  the  suspension. 

In  the  light  of  all  the  above,  we  consider  that  the  standard 
shipping  container  vibration  test  should  cover  only  the  frequency 
range  2-60  cycles  per  second. 
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CYCLING  vs.  RESONANCE  TESTING 


The  next  question  concerns  whether  the  test  should  be 
cycling,  resonance  or  combinations  thereof.  Practically  all 
tests  used  require  a  search  through  the  stipulated  frequency 
range  to  find  the  points  of  resonance  or  maximum  magnification. 
After  finding  resonances,  the  container  is  vibrated  at  these 
points  at  a  stipulated  amplitude  for  a  stipulated  period  of  time. 

The  major  exception  is  Picatinny  Arsenal  which  follows  the 
complex  schedule^  shown  in  Table  1. 

One  could  adopt  the  dangerous  premise  that  a  standard 
should  not  accept  the  Table  I  schedule  on  the  simple  grounds 
that  the  other  agencies  don^t  use  it.  Such  a  procedure  is,  of 
course,  scientific  nonsense  and  our  reasons  for  favoring  a 
simpler  approach  should  be  spelled  out.  We  favor  a  procedure 
which  involves  searching  out  the  fundamental  "resonant  frequency" 
and  vibrating  at  that  frequency  for  fairly  long  time  because; 

a.  We  are  evaluating  the  fatigue  resistance  of  the 

vibration  isolation  system. 

b*  Maximum  system  response  is  obtained  at  the 

fundamental  resonant  frequency.  Hence,  stress  levels  are  at 
the  maximum  which  can  be  obtained  with  motion  generating 
equipment  being  used. 

c.  The  cumulative  damage  hypotheses  indicate  that 
the  maximized  stresses  will,  for  the  shortest  test  time,  produce 
the  most  damage. 

d.  Reversing  c,  if  a  system  passes  the  test;  confidence 
in  service  reliability  of  the  system  is  obviously  maximized. 

e*.  While  it  is  true  that  one  can  select  the  wrong 
stress  or  time  parameter,  or  both,  further  studies  of  the 
environment  can  improve  the  validity  of  these  parameters.  Our 
theoretical  approach  supports  this  thesis. 


^The  schedule  reproduced  here  covers  only  ambient  conditions. 
Picatinny  also  requires  vibration  testing  with  different  time  constants 
at  temperature  extremes. 
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f.  All  other  things  being  equal,  the  simplest  approach 

is  to  be  preferred.  This  paper  indicates  that  a  simple  approach, 
is  feasible. 
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TEST  AMPLITUDES 


Having  disposed  of  the  questions  of  frequency  range  and 
type  of  test,  the  difficult  problems  of  controlling  amplitude  and 
time  to  simulate  the  effects  of  the  environment  now  follow.  For 
convenience,  amplitude  is  discussed  first  although  time  is  closely 
related.  Because  of  this  interrelation  we  will  discuss  amplitude 
as  a  target  optimum  and  discuss  possible  variations  in  amplitude 
and  time  later. 

Our  theoretical  study  indicates  that  there  is  no  input 
amplitude  function  (such  as  displacement-frequency  or,  accel¬ 
eration-frequency)  which  would  resiilt  in  equal  life  expenditure 
for  every  suspension  system.  Thus,  the  specification  writer's 
ideal  of  a  simple-test  applicable  to  all  suspensions  is  a  will-o'- 
the-wisp,  On  the  other  hand,  any  reduction  in  the  variety  of  tests 
currently  prescribed  is  an  obvious  desideratxxm. 

Equations  and  [44J  show  that  input  amplitude  can  be 

varied  depending  upon  the  response  characteristics  of  the  item 
on  its  suspension.  These  response  characteristics  (  On  )»  are,  in 
turn,  functions  of  the  acceleration  spectral  density,,5L  q. 
the  absence  of  statistical  measures  of  §4cu  have  a 

logical  basis  in  theory  for  stipulating  test  time  or  response  or  input 
amplitude.  While  it  should  be  possible  eventually  to  obtain  this 
information,  and  the  results  of  such  knowledge  acquisition  may  be 
totally  different  amplitudes  from  those  now  specified,  it  is  illogical 
to  force  vibration  testing  equipment  development  along  a  road,  the 
direction  and  length  of  which  cannot  begin  to  be  estimated  at  this 
time.  That  is  to  say,  we  have  no  right  to  ask  for  more  capacity 
than  is  available  until  we  know  we  need  it. 

In  furtherance  of  this  thought,  therefore,  we  investigated 
the  capacity  of  various  commercially  available  vibration  testers. 
Figures  11  through  16  show  the  capacity  of  various  machines  to 
move  the  several  weights  shown  through  various  amplitudes,*  One 
point  that  stands  out  above  all  others  is  that  there  is  only  one 
shaker  which  can  cover  all  the  amplitudes  and  frequencies  shown 


*The  drawings  are  not  cumulative  but,  rather,  represent  the 
maximum  capacity  of  the  existing  known  machines,  Obiyiously, 
large  machines  can  all  be  used  for  testing  small  packages. 


38 


in  Figure  2  to  be  currently  required  between  2-60  cps.  *  When 
amplitudes  are  restricted  to  those  producing  an  acceleration  of 
1.  3g  on  5000  lbs.  gross  weight,  the  Wyle  W-1000  type  can  per¬ 
form  at  all  frequencies.  **  Maximum  acceleration  output,  at 
5000  lbs.  gross  for  the  MB  HC-50  hydraulic  types  is  Ig, 

Obtaining  more  than  1  inch  double  amplitude  below  5 
cycles  per  second  is  a  very  expensive  endeavor.  As  it  now  stands, 
the  only  location  where  this  could  be  done  today  is  at  Norair  Division, 
Northrop  Corp.  ,  Hawthorne,  California***  which  would  add  con¬ 
siderably  to  the  cost  of  testing  a  container  designed  and  fabricated 
east  of  the  Rockies.  Lwet  us,  therefore,  look  a  little  more  closely 
at  the  requirements  of  those  specifications  which  ostensibly  call 
for  double  amplitudes  in  excess  of  1  inch.  Such  a  procedure  dis¬ 
closes  the  following: 

a.  MIL-P-9024A  and  MIL-E-4970A  show  the 
very  high  double  amplitudes  below  5  cps  identified 
in  Figure  2.  Both  MIL-P-9024A  and  MIL-E-4970 
permit  significant  deviation  as  exemplified  by  the 
following  from  paragraph  4.  6.  8.  1  of  MIL.-E-4970A: 

” . afa  suitable  sweep  frequency  rate  and 

with  a  vibratory  double  amplitude  or  a  peak  accel¬ 
eration  input  of  sufficient  magnitude  to  establish 
significant  resonant  frequencies.". 

b.  Picatinny  Arsenal' s  procedure  states  "When 
the  natural  frequency  lies  in  the  2-20  cps  range 
input  when  cycling  through  resonance  shall  be  reduced 
to  those  values  listed  in  the  resonance  test  for  2  cps 
below  to  2  cps  above  the  resonant  frequency".  Hence, 


*  This  is  the  Northrop  electro-hydraulic  system  of  which  the  Wyle 
company  is  a  licensee  for  the  shaker  components.  Northrop  has 
estimated  that  development  costs  for  their  system  were  around 
$700,  000.  00  and  that  cost  of  duplication  of  their  facility  would  be 
about  $300,  000.  00. 

**  An  informal  cost  estimate  from  Wyle  for  a  system  of  this  nature 
is  $48,  000  and  6  months  delivery. 

***  Northrop  has  expressed  willingness  to  make  time  available  on 
their  machine  to  others  and  has  stated  that,  once  accepted,  outside 
tasks  will  have  equal  priority  standing  with  in-house  work. 
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for  resonant  testing  in  the  low  range^  the  double 
amplitude  is  really  0.  13  inches. 

c.  Until  quite  recently,  there  were  no  commer¬ 

cially  available  equipments  which  could  shake  a 
container  through  an  approximately  8  inch  stroke. 

Hence,  most  previous  tests  have  been  run  at  1- 
inch  double  amplitude  obtainable  from  the  standard 
package  testers  covered  by  ASTM  D999-48T. 

In  view  of  the  foregoing,  practically  no  vibration  tests 
have  been  run  at  a  double  amplitude  in  excess  of  1.  0  inches  and 
yet  the  vast  number  of  tests  which  have  been  run  at  this  amplitude 
gives  an  indirect  measure  of  confidence  that  can  be  relied  upon  in 
a  very  large  number  of  cases.  *  We  have  machines  available  to  us 
which  will  vibrate  at  this  amplitude  between  2.  5  and  5  cps,  the 
L.  A.  B.  and  Gates  package  testers**.  With  minor  changes  in 
pulleys,  the  machines  can  be  modified  to  operate  between  2  and 
4,  5  cps.  As  will  be  brought  out  later,  4.  5  cps  (Ig  at  1  inch  double 
amplitude)  is  considered  high  enough.  Aside  from  the  fact  that 
most  tests  have  been  run  at  this  amplitude,  this  simple  change 
makes  it  possible  to  substitute  a  machine  costing  around  $5,  000 
max.  for  one  costing  almost  10  times  as  niuch.  This,  then, 
represents  our  recommendation  for  the  2  -  4.  5  cps  range. 

Over  4,  5  cps  one  must  choose  between  the  1.  5g  (stipulated 
by  Picatinny  up  to  7  cps)  the  1.  3g  of  MIL-P-21927  (NOrd),  MIL- 
P-9024A  and  MIL-E-4970  and  some  of  the  lesser  units.  Here  we 
take  the  view  that  the  maximum  practical  number  of  machines 
should  be  usable.  Looking  briefly  at  the  5000  lb,  capacity  curves 


*  There  is  no  way  one  can  avoid  the  implications  of  the  National 
Safe  Transit  Committee  program.  One  of  the  authors  has  personally 
conducted  such  vibration  tests  many  times.  Among  other  things  that 
have  been  predicted  successfully  and  cures  foxmd  in  the  laboratory 
are  vibration  damage  to  TV  sets  and  electronic  organs,  fatigue 
failures  of  an  overloaded  shock  mount  system,  fretting  corrosion 
damage  to  empty  aluminum  foil  dinner  trays,  improved  design  of 
hot  water  heaters,  causes  of  shipment  damage  to  large  electron 
tubes,  etc, ,  etc.  The  basic  NSTC  procedure  is  not  as  bad  as  some 
people  think.  When  properly  used  and  interpreted,  it  can  be  and  is 
a  powerful  tool, 

**  The  Gates  testers  are  not  identified  as  such  in  the  various  graphs. 
In  general,  what  L.  A.  B.  can  do,  Gates  can  do  and  at  a  competitive 
price. 
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(Figure  11)  we  find  an  awkward  gap  between  4.  5  cps  and  5  cps  at 
the  1.  3g  level  (assuming  that  the  LiAB  machine  has  been  modified 
to  get  down  to  2  cps)  where  only  one  machine  in  existence  can 
produce  the  necessary  frequency  and  amplitudes.  This  same  gap 
exists  all  the  way  down  to  100  pound  capacity  machines.  *  Over 
5  cps,  at  5000  lbs.  MB  CTO  and  Ling  246  have  the  theoretical 
capacity.  If  we  drop  the  input  to  Ig,  however,  we  bring  into  the 
picture  the  MB  HC50  family  of  hydraulic  shakers  plus  Royal  Jet's 
special  and,  above  5  cps,  the  whole  family  of  MB  electromechanical 
shakers.  In  other  words,  stipulating  Ig  above  4.  5  or  5  cycles  per 
second  results  in  maximum  machine  availability  for  performing 
the  necessary  tests.  Of  particular  importance,  the  MB  line  of 
machines  is  made  usable.  Regardless  of  any  merits  or  demerits 
of  these  particular  devices,  the  MB  concern  has  sold  more  of 
their  machines  than  any  other  company,  partly  because  they  were 
among  the  very  first  in  the  field.  In  the  absence  of  a  statistically 
sound  description  of  the  environment,  we  can  see  no  compelling 
reason  for  eliminating  the  majority  of  available  devices  from 
consideration  just  because  of  an  arbitrary  interpretation  of  that 
environment.  We  therefore,  recommend  an  equally  arbitrary 
decision  setting  the  amplitude  at  Ig  between  4,  5  and  60  cps.  The 
complete  frequency  amplitude  spectrum  recommended  as  a  possible 
standard  is  shown  in  Figure  17. 


*  This  gap  is  not  always  present.  Thus,  it  is  possible  to  modify  a 
"package  fester"  so  that  it  will  cover  both  2  and  5  cps.  Further  at 
5000  lbs.  the  Ling  and  MB  machines  might  be  run  down  that  extra 
1/2  cycle.  If  the  load  is  considerably  less  than  maximum  capacity, 
you  do  have  a  chance  of  covering  the  gap  with  a  machine  of  consider¬ 
ably  greater  capacity  than  would  normally  be  rated  for  the  job. 
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TEST  TIMES 


Before  discussing  possible  equivalent  variances  in  annpli- 
tudes  to  further  broaden  the  range,  let  us  first  examine  the  time 
parameter.  Our  theoretical  studies  disclose  no  basis,  with 
currently  available  data,  for  selecting  an  optimum  test  time  for 
the  amplitudes  selected.  We  proceed,  therefore,  to  reasoning 
by  analogy. 

Most  of  the  tests  require  1/Z  hour  at  resonance  but  at  the 
1.  3  g  intensity  level.  In  the  range  of  frequencies  covered  by  the 
National  Safe  Transit  Committee  procedures,  and  with  which  the 
author  has  had  the  most  experience,  it  is  considered  that  a  sound 
system  should  be  capable  of  withstanding’ 1  hour  at  "resonance". 
Commercial  tests  are  based  on  this  time  scale  at  a  1"  double 
amplitude  and  it  hardly  seems  reasonable  to  require  less  in  guided 
missile  containers.  The  authors  own  personal  experience  is  that 
reliable  systems  will  last  almost  indefinitely  at  these  amplitudes 
and  frequencies  whereas  many  which  are  marginal  require  more 
than  30  minutes  to  fail*. 

Inasmuch  as  we  have  recommended  a  Ig  intensity  in  lieu 
of  the  more  commonly  specified  1.  3g  at  frequencies  in  excess  of 
4,  5  cycles  per  second  some  increase  in  the  time  at  resonance 
over  the  30  minutfes  usually  stipulated  for  this  range  appears 
warranted.  For  conformity  with  the  low  frequency  test  require¬ 
ment,  1  hour  is  recommended**. 


*  A  specific  case  of  a  hi-fi  unit  is  recalled  in  which  the  chassis 
was  supported  on  soft  rubber  mounts.  An  alarming  number  of 
instances  of  transformer  mounting  lug  failure  had  been  reported. 

In  the  laboratory  after  about  40  minutes  at  3.5  cycles  1"  double 
amplitude,  fatigue  failure  of  the  lugs  occurred.  The  transformer 
then  "tooh  off"  and  "cleaned  out"  the  inside  of  the  hi-fi  set.  Appro¬ 
priate  corrections  were  then  devised  and  no  failures  occurred  in 
1  hr.  No  further  failures  in  shipment  were  then  encountered.  While 
this  was  actually  a  case  of  faulty  unit  design,  the  principle  disclosed 
is  valid  for  the  package, 

**  Since  we  are  also  recommending  vibrating  only  at  the  fundamental 
resonance,  this  increase  in  time  at  this  point  will  still  result  in  a 
net  decrease  in  overall  test  time. 
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Heretofore,  we  have  been  discussing  testing  as  though  the 
entire  container  could  be  moved  in  conformity  with  the  input 
motions  stipulated.  There  are,  however,  many  containers  that 
simply  will  not  fit  on  existing  equipment  without  substantial, 
and  quite  often  mechanically  3ignificant,  modifications  in  the 
tester.  Figure  18  summarizes  some  of  the  problems  associated 
with  testing  when  machines  of  adequate  size  and  force  capacity 
are  not  available.  If  a  practical  test  procedure  is  to  be  written, 
then  the  practical  test  procedure  must  recognize  the  practical 
difficulties  of  compliance.  At  the  same  time,  it  is  felt  that  a 
practical  test  should  tahe  into  account  the  possibility  of  further 
increasing  severity  and  reducing  test  time.  As  pointed  out 
earlier,  current  specifications  permit  reduction  in  amplitude  but 
do  not  make  any  attempt  to  modify  the  time  parameter.  Equations 
[43j  and  [44]  ,  on  the  other  hand,  indicate  that  change  in  test 
time  with  change  in  amplitude  is  a  logical  necessity.  In  the 
absence,  however,  of  appropriate  data  concerning  proper  values 
for  the  slope  of  the  S-N  curve  (3),  direct  use  of  these  equations 
is  not  possible  at  this  time. 

If  we  are  to  propose  a  standard  at  this  time  which  takes 
into  account  the  practical  difficulties  of  complying  with  a 
stipulated  amplitude,  we  must  again  be  arbitrary  while  still 
remaining  consistent  with  our  theoretical  study.  In  any  case, 
a  specification  must  be  consistent  with  itself  even  though  it  may 
not  be  consistent  with  reality.  One  way  in  which  time  might 
reasonably  be  varied  is  to  equate  the  .total  input. work.  Since 
work  for  a  given  sine  wave  cycle*is  conveniently  measured  by 
the  mean  square  amplitude,  we  find  that  total  work  is  the  same 
when  the  ratio  of  test  times  is  inversely  proportional  to  the  mean 
square  amplitudes,  or 

t'r  , 

tv  VX'CtJ 


This  expression  is  the  same  as  the  similar  expression  for 
the  log-log  S-N  curve,  simply  through  substitution  of  3  =  1/2 
in  equation  0^  .  Thus,  our  criterion  of  maintaining  some 

consistency  with  the  theoretical  approach  is  satisfied. 

^  By  the  cumulative  damage  hypotheses  the  relation  given 
above  can  only  be  valid  in  the  range  X  least  large  enbugh 

to  just  exceed  the  infinite  number  of  stresses  line  on  the  S-N  curve 
to  the  maximum  one  time  stress  the  system  is  capable  of 


TEST  FIXTURE 


FIGURE  18 

SOME  VARIATIONS  ON  BRUTE  FORCE  VIBRATION  TESTING 


withstanding.  Since  these  values  are  not  known  these  limits  are 
somewhat  theoretical.  We  take  refuge  in  analysis  of  Figure  18 
which  tells  us,  in  effect,  that  it  is  not  essential  that  the  input 
amplitude  be  less  than  1/2  the  stipulated  value*.  For  a  high 
limit  we  can  take  refuge  in  the  opinion  expressed  by  one  industry 
responder  that  10  minutes  at  resonance  with  a  high  amplitude 
input  should  normally  be  sufficient.  A  plot  of  the  relation  given 
above  in  terms  of  tv  60  npinutes  and  within  the  limitsXy  X 

tff  > 

.j-  s  10  is  shown  in  Figure  19  which  is  recommended  for 
adoption. 


*  In  fairness  it  is  not  proven  that  this  is  a  valid  lower  limit.  The 
author  has  been  ’’burned'*  by  a  system  vibrated  at  approximately 
3  cps  with  one  end  on  an  LAB  package  tester  as  in  scheme  D  of  Figure 
18.  With  this  low  natural  frequency  the  shear  mount  system  had 
auxiliary  friction  dampers  such  that  transmissibility  was  less  than 
3,  After  considerable  cut  and  try,  including  use  of  teflon  inserts 
in  the  damper  members,  it  was  possible  to  get  the  system  to  pass 
the  contract  requirements  of  1/2  hour  input  at  each  end.  When  the 
system  was  plat:ed  in  serv  ice,  however,  the  dampers  heated  up, 
seized  and  then  parted  company.  Many  instances  of  broken  mounts 
were  found.  Eventually,  the  system  had  to  be  completely  redesigned 
to  a  different  natural  frequency  before  a  serviceable  container  was 
achieved.  In  the  questionnaire  results  many  of  the  responders 
regarded  this  particular  situation  as  a  freak.  But  what  good  is  a 
specification  if  it  doesn’t  at  least  screen  out  the  freaks?  All  we 
know  for  sure  is  that  1  hour  at  half  the  stipulated  amplitude  wasn’t 
enough  and  that  a  four  hour  test  might  have  caught  it.  Or,  four 
hour  testing  certainly  should  have  given  us  a  better  chance  to  find 
this  element  of  unreliability  than  the  1  hour  test  which  failed  to 
prove  it. 
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FIGURE  19 

SUGGESTED  TIME  OF 
TEST  AT  AMPLITUDES 
DIFFERING  FROM  THOSE 
GIVEN  IN  FIGURE  17 
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DIRECTION  OF  APPLICATION 


Practically  all  of  the  test  procedures  analyzed  require  that 
input  vibrations  be  applied  to  the  three  orthogonal  axes  of  the 
container  on  the  assumption,  apparently,  that  the  unit  has  equal 
chance  of  being  stowed  in  the  carrier  on  any  face.  The  mlhute, 
however,  that  the  container  has  so  much  as  a  ^‘This  Side  Up^* 
marked  on  it  the  chances  are  no  longer  equal  that  any  specific 
container  will  be  stored  on  any  face*.  The  moment  a  container 
has  a  definite  base  defined  in  some  manner,  such  as  by  skids, 
etc,  ,  carriage  on  any  other  face  becomes  a  difficult  achievement. 

From  the  standpoint  of  the  designer ^  the  possibility  of 
vibration  on  any  face  requires  careful  design  analysis,  frequently 
imposing  complexities  in  the  design  vhich  might  otherwise  be  un¬ 
necessary,  Almost  all  knowledgeable  designers  oi  large  containers, 
therefore,  seek  deviation  from  the  requirement  in  one  fashion  or 
another  and,  in  most  cases,  this  deviation  is  granted,  tacitly,  if 
not  explicitly.  Since  the  deviation  appears  to  have  been  granted 
more  often  than  not,  correct  specification  practice  dictates  placing 
all  designers  on  equal  footing.  Stipulation  that  vibration  be  applied 
only  through  the  base  on  skidded  containers  appears  to  be  an  obvious 
requirement,  ^ 


*In  spite  of  our  usual  cynicism  in  these  matters,  some  people  do 
read  and  heed.  Incidentally,  a  carrier  can  be  held  liable  for  damaged 
goods  if  it  can  be  shown  he  ignored  such'  a  warning. 
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FASTENING  TO  THE  VIBRATION  TABLE 


There  is  considerable  varia.Uon  in  stipulating  the  method 
of  fastening  the  container  to  the  vibration  sou’ ce,  MIL-W-21927 
(NOrd)  prescribes  tying  down  as  it  would  be  in  the  rail  car  while 
other  specifications  tend  towards  use  of  the  jfarase  "rigidly 
attached".  While  ideal  infinite  rigidity  is  impossible,  vibration 
test  technicians  much  prefer  the  prettier  results  obtained  from 
a  "rigid"  attachment.  This  occasionally  produces  the  odd  result 
that  the  tie  downs  are  more  massive  than  the  container  I  In  many 
cases,  the  tie  downs  are  considerably  more  expensive  than  what 
would  be  used  in  practice.  This  element  of  unreality  is  further 
emphasized  by  a  quick  glance  at  most  carloading  practices  for 
heavy  containers.  In  most  blocking  arrangements,  the  restraints 
are  primarily  fore  and  aft  with  vertical  restraints  consisting 
mostly  of  far  from  rigid  steel  strapping.  For  smaller  containers 
the  vertical  restraint  usually  consists  only  of  the  weight  of  super¬ 
imposed  containers  plus  friction  of  adjacent  containers.  Possibility 
always  exists  that  some  significant  portion  of  the  journey  could 
well  be  without  any  restraint. 

A  completely  unrestrained  container,  on  the  other  hand, 
could  "take  off".  Aside  from  the  fact  that  falling  off  the  vibration 
rig  could  be  disastrous,  the  element  of  personnel  safety  enters. 

Now  the  NSTC  test  procedure  does  not  actually  call  for 
testing  at  the  point  at  which  maximum  resonance  occurs.  The 
procedure  stipulates  that  one  increases  the  frequency  until  such 
time  as  the  resultant  causes  the  outer  container  to  lift  off  the 
table  approximately  1/32  inch.  This  may  or  may  not  be  the  frequency 
producing  maximum  transmissibility  but  it  is  not  far  removed  from 
it.  From  the  number  of  tests  conducted,  this  slight  extra  input  at 
each  cycle  seems  to  improve  the  confidence  in  the  test  results.  In 
any  case,  the  slight  additional  motion,  while  defying  generalized 
analysis,  can  be  considered  as  a  possible  allowable  increased 
motion  which  does  not  "ruin"  the  test  even  though  it  certainly 
will  confuse  the  Record,  We,  therefore,  make  the  suggestion  that 
this  amount  of  motion  be  allowed.  It  should  be  noted  in  passing 
that  holding  acceleration  to  1  g  in  lieu  of  1.  3  g,  greatly  reduces 
the  strength  requirements  for  hold  downs  and,  therefore,  contributes 
to  a  lesser  cost  of  testing. 
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PROPOSED  LANGUAGE  FOR  A  STANDARD  TEST 


With  the  above  reasoning  in  m^ndj,  it  is  possible  to  write  a 
test  procedure  requirement  which  could  be  adopted  as  standard 
and  which  approaches  our  self  imposed  criteria.  This  proposed 
language  follows, 

Basic  Test 

The  container,  loaded  with  its  contents  or  with  a 
dummy  of  equal  weight  and  equivalent  principal  moments  of 
inertia,  shall  be  vibrated  in  such  a  fashion  as  to  discover  the 
frequency  producing  the  maximum  transmissibility  within  the 
frequency  range  of  2  to  60  cycles  per  second.  The  container 
shall  then  be  suitably  fastened  {maximum  permissible  motion 
freedom  approximately  l/32nd  of  an  inch)  to  the  table  of  a 
vibration  testing  device  producing  an  essentially  sine  wave 
motion  and  vibrated  for  one  hour  at  the  frequency  producing 
’maximum  transmissibility  and  at  the  double  amplitude  shown 
in  Figure  17. 


The  test  shall  be  applied  to  each  of  the  three 
principal  orthogonal  axes  of  the  shipping  container  except 
that  where  containers  are  equipped  with  a  definite  base,  defined 
by  skids  or  equivalent,  vibration  inpu^^s  shall  be  applied  only 
through  such  base.  This  test  shall  preferably  be  applied  just 
before  any  required  shock  tests. 

Permissible  Alternates 


a.  When  a  device  such  as  one  meeting  ASTM 
specification  D996  is  used  for  frequencies  less  than  5 
cycles  per  second,  the  horizontal  motion  inherent  in 
such  machines  is  acceptable.  Wherever  possible,  how¬ 
ever,  machines  shall  be  adjusted  to  produce  ’’vertical- 
linear"  motion. 

b.  When  testing  devices  are  not  available  to 
produce  the  required  double  amplitude  over  the  entire 
container,  alternate  methods  will  be  acceptable  provided 
the  motion  at  the  locations  of  mounts  (or  the  centroids 

of  each  cushion  on  either  side  of  the  item  center  of  gravity) 
is  the  same  as  that  specified.  A  test  procedure  which 
places  one  end  of  the  container  on  a  pivot  while  the  other  end 
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is  moved  is  also  acceptable  provided  the  position  of 
the  vibration  generator  is  changed  from  one  end  of 
the  container  to  the  other  as  rapidly  as  possible  half¬ 
way  through  the  test. 

c.  Variations  in  stipulated  double  amplitudes 

are  permitted  provided  that  the  time  of  teat  be  adjusted 
in  accordance  with  Figure  19.  The  minimum  time  at 
resonance  shall  be  10  mihutes  and  the  maximum  time  at 
resonance  shall  be  4  hours. 
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DISCUSSION 


It  is  regretted  that  a  search  for  transportation  vibration 
data  in  a  form  useful  for  computing  in  accordance  with  the  theories 
advanced  herein  produced  negative  results.  Most  literature  data  are 
already  reduced  to  a  form  useful  for  the  particular  study  in  hand. 

If  original  records  are  still  available,  it  might  be  possible  to  reduce 
the  original  data  in  the  form  needed  here.  It  is  considered  of 
fundamental  importance  that; 

a.  Every  effort  be  made  to  salvage  as  much  random 
data  from  still  existing  records  as  possible. 

b.  Future  records  be  analyzed  for  whatever  purposes 
desired  but  also  be  analyzed  for  contribution  to  a 
proper  statistical  description  of  the  environment. 

Achieving  a  useful  description  of  the  environment  is  of 
fundamental  importance  regardless  of  the  ultimate  fate  of  the  theory 
advanced  herein.  The  theory  cannot  be  completely  discounted  until 
such  time  as  the  environment  is  shown  not  to  conform  to  our 
assumptions. 

We  have  also,  however,  developed  a  theory  of  possible 
behavior  based  on  the  presumption  that  cumulative  fatigue  hypotheses 
apply  to  package  shock  and  vibration  isolation  devices  and  materials. 

This  assumption  lacks  any  experimental  verification  except  for  steel 
springs.  The  scope  of  this  study  precluded  such  experimentation 
which,  in  any  case,  should  not  have  been  undertaken  until  a  theory 
needing  verification  existed.  Now  that  the  theory  exists  it  is  considered 
a  most  urgent  necessity  that  work  be  undertaken  to  verify,  on  a  spot 
check  basis  as  a  minimum,  whether  or  not  shock  isolation  systems 
obey  cumulative  damage  hypotheses. 

Complete  verification  of  all  facets  of  our  theory  cannot  come 
immediately.  After  studying  the  various  test  methods  currently  specified 
the  only  reasonable  conclusion  is  that  some  form  of  standardization  is 
an  immediate  requirement.  We  have,  therefore,  sought  a  logical  means 
of  cutting  through  the  fog  to  some  basis  which,  although  of  necessity 
rather  arbitrary,  recognizes  the  practical  limitations  of  conducting  any 
such  test  and  is  consistent  with  at  least  one  of  the  theoretical  lines 
developed  herein.  To  the  best  of  our  knowledge,  no  other  testing  method 
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proposal  can  make  such  a  claim  while  $till  maintaining  the  virtue 
of  compliance  ability  on  a  broad  national  basis.  Accordingly, 
it  is  urged  that  the  test  method  proposed  be  adopted  as  an  interim 
measure.  The  fact  that  this  interim  may  last  a  number  of  years 
simply  buttresses  any  arguments  against  waiting  for  a  perfected 
method. 
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APPENDIX  I 


MATHEMATICAL  DERIVATIONS 


In  the  interest  of  maintaining  continuity  of  reasoning  some 
rather  tedious  mathematical  operations  have  been  omitted  from 
the  main  body  of  the  report.  This  appendix  details  these  operations 
for  the  record. 


Derivation  of  Equation  Bzl  ■ 
Equation  Da  reads 

while  equation  m  reads 

Sia,to5)  =  B t  / 

By  direct  substitution  (A-1)  becomes 

a 

^  t  ULJ  ^ 

a; 


_ Qi  ^  SiJ(o\i 


(A-l) 


(A -2) 


O  (A- 3) 


Equation  £7'2  gives  a  value  of 


CO 


I  to 


ri 


Squaring  and  substituting  into  (A- 3)  leads  to  the 
expression 


(A -4) 


rational 
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For  simplicity  then  cjcOsCO^cl  V  and  apply 

CO'  ^ 


the  further  reasonable  restrictions  that  ^0  i8i<^  2  and  .  . 

^  .  For  a  given  set  of  conditionsto^^CtuS-nd  Suouff^^n) 

will  be  constants.  Hence  (A-5)  goes  over  to 


(A-6) 


be  seen  that  we  have  a  rational  function  of  ^  which  can  be 
integrated  if  expressed  in  terms  of  partial  fractions. 

The  left  hand  member  of  the  denominator  in  the  integrand 
of  equation  (A-6)  expands  and  then  factors  into 


{A-7) 


where 


Similarly  the  right  hand  term  of  the  s  ame  denominator 
expands  and  factors  into 

(\-  et  X  =  X\x"-(  bt-  1  ,A.8, 


where 


»(x  i/Jx  +  i)^x^/3x + 


Hence 


y3- 

f,(x^  will  assume  the  form 


A  +  C  -K  P.  4-Ad>N  X  ■  ^  P+Q  X — 

”  X'V/3)(+!  x%^x+l 
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Now,  let 

E=  KtD  *C  +=r  (i-  k),£  =  D  (C-4)^ 

E^=  A  +C ,  G,=N  tQ,  G^=  M  4  Pt  a  (Q  -  N) , 

G=Q+/V+/?(P-/M)  and  sD  + 

then  multiplying  the  left  hand  sections  (in  braces)  of  the  numerators 
of  the  two  fractions  and  substituting,  we  obtain  ‘ 


Expanding  and  equating  to  the  appropriate  coefficients  of  yC  we 
obtain 
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A- 12  (a),  (b),  (g)  and  (h)  give  us 

Ei  +  G|.0,  E,= 

E  a+  "  ^a. 

E 1+  63  =  Q  ^3*  ~ 

Substituting  in  A- 12  (c),  we  have 

t.(6?-4r)=0 
E  t  -  G(=  o 

A-12  (d)  reduces  to 

E.-G.- 

A-12  (e)  gives 

^3=  ^3=0 

and  A- 12  (f)  gives 


'S 


&1-4  r 


E.=-g 


I 


^  Bt-4r 

Dropping  the  intermediate  substitutions  we  write 

+  D  =  0 

E^.AtC+.^rfb-KV  ^ 


E=A-.-C=  — I - , 

which  leads  directly  to 


Bi-4T 


D=-K-  4~C^-  I 

also,  DtK+o^(/4-C)=0 

4-Ca  _ i- _ 

a(Bt-4T")  ■ 


(A- 13) 


(A- 14) 

{A- 15) 


{A-16) 


{A-17) 


(A-18) 


(A- 19) 


(A-20) 
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Similarly,  we  find 

N=  -Q  = 
P=M  - 


Using  these  values  of  the  arbitrary 
becomes  r 


aitrary  constai 


(A-21) 


nts  equlation  (A -9) 


_J _  uV^  li-^X 

lx+o£')<4i  xS/?x+i 

^  (A-ZZ) 


Let  .  Then  the  integrand  becomes 


'  Xtt^x+l  y  y^-/x4i 

Integrating  each  part  using  Pierce's  formulas  and  [72]  ,  we  find 


-l±lSl2L_clx«^/- _ ^ _  I  rT\  j  xdx 

X^+o^x+l  ^ 


T 

=LW' 


<5-  C5C  ^ 


.fi*|(xfe(i+j)-|CW 


(A-24) 
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Similarly 


r 


^  ('xVSs^ti) 


(A-26) 


(A-27) 


Hence 


f 


f,  Cx)clx  St 


m 

Ed: 

m 


X^t£t  I 
X^-oC)H-i 


-.m 

20 


X!±£)(i! . 


2(BtUn 


IjJqA 

jj 


-  -g  /ii£.^  j±iS!\ 

'  8(B^4nV  / 

4B.?fB+2y) 

& 


Therefore 


O  2ir  U  V.  ni  4b.^5(6^+£J) 


(A-29) 


Which  is  the  form  expressed  for  equation  [22]  in  the  main 
body  of  the  report. 

B.  Derivation  of  Equations  C35]  and  [36]  . 

Equation  (^31^  reads 

(B-i) 


while  [32]  and  [34]  are 


“''“■sV'ctT®  ^■'=‘'‘2 


(B-2) 


h  Xp 

,  '-ET  -r 

IT*®  € 

Hu 


Direct  substitution  produces 


p/4.^  -  ® 

ECt)-  ‘^j 


r 


Xp  Xp 
B  “  Wr,  I 

Xpe  dxp 


te. 


Let 


I 

h 


B  X 


2a 


from  which  we  immediately  obtain 


dXp  =  jBOnd^ 

^pdXp  —  dOn.'^c}^ 


J 


(B-3) 


(B-4) 


(B-5) 
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Let  us  also  introduce  an  appropriate  value  for  the  endurance 
limit  motion  X ,  such  thatXe  — /Sl'O'n  K  •  becomes 


-A  r^Tt 


Ir 


(B-6) 


Now  let  and  we  write  the  integrand  a 

b 


(B-7) 


In  evaluating  ,  complete  the  square 


For  simplicity,  let 
ci  =  ci  ij.Then 


(B-8) 


^  so  thafV-U +-3»- and 


TT 


Sl 


CL 

c 


-vjSL 


(B-9) 
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J*3  ^ 


-1 


— *o 


Pierce,  492: 

T  =s^ir*^/^ 

-1-4  a.  3l  4  / 

Substituting  in  [B-,]  we  have 


(B-10) 


(B-11) 


I- 


(B-12) 


(B-13) 


f  V  ^  J 

Turning  to  J.  g^we  again  complete  the  square 

which,  by  making  the  same  substitutions  as  before  becomes 

:-4" /2-^d  u , 


I  -  c 


A. 

Therefore 


-tit 

A\ 


I..., 


5 

(B-15) 


^  -o 

Multiplying  X(,^y  transposes  the  integrand  of  X^ 

into  the  probability  integral 


Jience, 


1  =  e 

•*•3, 


4- 


ja^. 


-K*'  r” 

I  -e  + 

X 


(B-17) 


Values 
page 


are  tabulated  in  the  literature,  e.  g.  Pierce 


14) 
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Substituting  in  we  have,  getting  rid  of  the  artificial 

constant  Q,  ,  ; 

a.  Sk 

"  'U 


r 


3- 


ag; 


l-e  + 

*—  i/H 


ctfr 


C*|(K) 


Ort. 


6 

C“T  — '  ^ 

B-6J  can  be  written  as 


18U 

X 


-1C 


t 

which  is  the  form  of  equation 


(B-19) 


It  can  also  readily  bo  seen  that  whenlC  is  very  small, 
or  non-existent,  [b-17]  ,  will  simplify  to 


■%5r" 


JL 

SI 


(B^20) 


which  is  the  form  given  for  equation  r36j 
C.  Derivation  of  Equations  [4^  and  [421  . 

The  equation  for  a  log-log  S-N  relation  is 

jLecj  Xpw  ^4?^  A-B  X'p^X^(C-l) 

’k 


From  which  hty^  = 


_A 

X 


P 


or 


/pAi 

A  ) 


(C-2) 
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Substituting  in  equation  we  have 


ECtI-  -%V 

aTT^  On 


(C-  3) 


'He 


Introducing  the  notation^^^-  and  its  corollaries  (C-3)  becomes 


B(ty- 


Stta^ 


-  ^nt 
TT 


(^)-k J. 


I 


r  w = £ 


a 

.  --■»  d, 

-^O 


a. 


Now,  the  complete  Gamma  function  in  ^  is  defined  by  the  relations 


and 


while  the  incomplete  Gamma  function  is  defined  by 

X(n,K)*2ry'’''e'^cj^ 


(C-5) 


(C-6) 


It  is  evident  that  the  integrand  of  equation  (C-4)  can  be  expressed 
in  terms  of  Gamma  fxinctions  with  rkai'_4  J  . 


SB 


Then,  we  have 


7 


i+i  -3® 


4-  r 


'+4 

e 


rO'*’ 20) "  (i ■'■^1  k)j 


.  niMi 

2 


1  ^6+A^* 

re -4) 


(C-7) 
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from  which  (  C-4  )  becomes 


l_  Y('\+-h'  K\ 


(C-8) 


The  life  expenditure  function  for^  sinusoidal  vibration  test 
with  a  log-log  S-N  curve  is  written  rather  simply  as 


E(t> 


-LA 

lx 


SlTT  V  . 

Combining  and  solving  for  time,  we  have 


\ -L 


tr  . 


=  Xn7fe(7r- 


ae 


(C-9) 


(C-10) 


which  is  the  form  given  for  equation  m- 

It  is  directly  evident  that,  if  K  very^mall^r  noij- existent, 


one  can  ignore  the  contribution  of  the  integral 
the  relation  of  (C-10)  can  be  simplified  to  ®  w 

C  T  -  '  ' 


B  e 


r 

(C-ll) 


and 


which  is  the  form  given  for  equation  [42]  . 


To  compute  solutions,  it  is  only  necessary  to  note  that  tables 
and  curves  of  (  )i^'^and  related  functions  can  be  found  in  Jahnke- 
Emde  who  also  refer  to  other  extensive  tabulations  of  these  functions. 
Jahnke-Emde  also  give  curves  of  the  ratio  If  (  H  1  1  V  *  while 

r  (h  +  h 

Pearson  has  compiled  extensive  tables  of  the  incomplete  Gamma 
function.  Hence  our  solutions  are  in  computable  form. 
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APPENDIX  II 


Glossary  of  Mathematical  Notation 


A 


empirical  constant  in  e^^ion  [34]  and  used  in 
derivation  of  equation 


defined  by  equation 
Fourier  transform 


w 

of  > 


^at 

Appendix  lA. 


is  the 


test  input  acceleration  amplitude  corresponding  to 
sinusodial  test  input  motion  X 

^  the  acceleration  of  motion 

c)t=- 

absolute  response  acceleration  corresponding  to 
^  '  absolute  response  motion  ^  o  (  "b  ^ 

Q,^  ^'^'^input  acceleration 


a 


-S. 


mean  square  absolute  acceleration  response 


Oh  used  in  |^-7^  ^ 


the  constant  first  used  in  equating  [34j  ,  which  is 

the  motion  constant  equivalent  to  the  .slope  of  the 
S-N  curve. 


C. 

Cc 

c 

D 

A 


the  bandwidth*G^f.^  0)*(^  is  the  frequency  difference 
between  the  points  on  each  side  of  the  peak  at  which 
the  ordinate  of  the  curve  is  half  the  peak  value. 

damping  coefficient,  lbs.  sec.  /in. 

critical  damping  coefficient. 

Miner's  constant,  ratio  of  stress  cycles  occurring  to 
stress  cycles  to  failure  (ostensibly  1,*0).  Also  used 
in  derivation  of  equation  ^22^  Appendix  lA. 

used  in  derivation  of  equation  Appendix  lA. 

experimentally  determined  exponent  in  the  Corten 
and  Dolan  non-linear  fatigue  constant. 
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=  the  error  function  of  K  sometimes 
"obability  integral. 


b.  Life  expenditure  in  Miner's  hypothesis 
^'^'^life  expenditure  as  a  function  of  time. 

1^^  ^'^'y.ife  expenditure  time  function  in  sinusoidal  test 
E.  I )  E  derivation  of  equation  Appendix  lA. 

"P  frequency,  cycles  per  second, 

P ^  natural  frequency,  cycles  per  second. 

^X^used  in  appendix  lA. 

Q  number  of  multiples  of  acceleration  of  gravity, 

jG-j  derivation  of  equation  Appendix  lA, 

M  complex  frequency  respon^  or  magnification 
factor  defined  by  equation  I  Sj  . 

H  height  of  drop  equation  [46]  ,  [47] 

g  , J- ^  symbols  to  identify  integrals  in  Appendix  1 

Kv  used  in  derivation  of  equation  [2Z-}  Appendix  lA. 

spring  rate  of  suspension,  also  in 

equation  |b-6*]  ^ 

ied  in  derivation  of  equation  [zz]  Appendix  lA. 


In 


use 


mass 


fs/  used  in  derivation  of  equation  [zz]  Appendix  lA. 
fatigue  life  at  maximum  stress  (Corten  and  Dolan) 
number  of  cycles  to  failure  at  amplitude  y  p  . 
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fatigue  life  at  stress  levelh  used  in  Miner's 
hypothesis  characterizing  the  linear  damage 
accumulation  in  fatigue  considerations. 

used  in  Corten  and  Dolan  s  non-linear  damage 
accumulation  proposal  as  the  total  number  of 
cycles  to  failure. 


mJ,  acutal  number  of  cycles  at  a  given  stress  level 
in  Miner's  hypothesis 

P  used  in  derivation  of  equation  Appendix  lA. 

p^probability  density  ofX  p  defined  by  equation 
Q  used  in  derivation  of  equation  Pzzl  Appendix  lA, 
Q  sometimes  called  quality  factor  0^ - where 


is  damping  ratio < 

autocorrelation  function  defined  in  a  footnote  to 
equation  M 

endurance  limit,  i.  e.,  stress  at  which  fatigue  is 
infinite 


stress  level 

Ipower  spectral  density  shows  the  power  content 
of  the  vibratory  components  as  a  function  of 
their  circular  frequency.  CO, 


power  spectral  density  of  the  input  motion. 
vDox  (^CO  J  power  spectral  density  of  absolute  response, 
Sq  (CO'^  power  spectral  density  of  relative  response. 

acceleration  spectral  density. 

So  Q,  spectral  density  of  the  absolute 


response. 
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f 


n'  input  acceleration  spectral  density  corresponding 
to  the  natural  frequency  of  the  suspension 


a  particular  value  of  time 

the  transmissibility  ratio  between  complex 
functions,  defined  by  equation  {7j 

time  coordinate 

test  time 

test  time  corresponding  to  response  amplitude 
test  time  corresponding  to  response  amplitude 


^nT 


w 

X 


nr  > 


weight 


see  equations 


amplitude  (distance  from  rest  position  to  peak)  inches 


mean  square  amplitude  defined  by  equation 


mean  square  relative  amplitude 


function  representing  motion  in  which  X'  denotes 
displacement  and  "tj  indicates  X  is  a  function 
of  time. 


X 


p  peak  value  of  X  ,  or  an  ordinate  of  the  envelope 
curve. 


X i^pu^t  motion  coordinate 

Xq  (""t  Jab  solute  response  motion  coordinate 

X,(t)  relative  response  motion  coordinate 

I X  complex  functions  of  the  circular  frequency 

see  equation  [0 

amplitude  of  a  sinusoidal  test  motion 
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Uniform  response  amplitude  to  ainusodial  teat 
motion. 

test  input  amplitudes 

test  relative  response  amplitudes,  see  equations  [4^ 

&l] 


Xr  amplitude  of  sinusodial  test  motion 

^  mean  square  relative  response  amplitude  in  test, 
amplitude  corresponding  to  endurance  limit. 


oC 

oT 


oC-. 


/3 


used  in  evaluating  integral  in  Appendix  I, 

ZX  Xp 

used  in  Appendix  A-? 

used  as  a  constant  in  semi-logarithmic  stress 
function. 


s 


percentages  of  cycles  at  stresses  v_)  ^  (Corten 
and  Dolan). 


used  as  a  constant  in  semi-logarithmic  stress 
function. 


C, 


/3  used  in  Appendix  A-8. 

complete  Gamma  function, 
the  incomplete  Gamma  function, 
ratio  of  damping  to  critical  damping  ^  ~  Q  ^ 
a  called  root  mean  square  and  also  standard  deviation. 

f  ^ 

O  r.m^s.  value  of  response  motion  ( CX  »  X 

^  I  ^  ' 

CO  circular  frequency,  radians  per  second 
natural  frequency,  radians  per  second 
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